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(54) Adaptive antenna 

(57) An adaptive antenna arrangement including a plurality of antenna elements in an array, a beam-forming 
network (20), one of the antenna elements providing a primary signal y(t), the remaining elements 
providing auxiliary signal inputs x(t) to the beam-forming network, a beam pattern controller (21) to which 
the auxiliary signals are applied, the controller being adapted to form amplitude and phase weights to be 
supplied to the beam-forming network whereby the array beam pattern is continuously adjusted to contain 
nulls which track the bearings of unwanted received signals, and means for combining the output of the 
beam-forming network with the primary signal, characterised in that the controller performs a so-called 
Q-R decomposition of the auxiliary signals to determine the weights to be applied to the beam-forming 
network and the array residual signal e(t) is fed back to the controller to modify the weights to correct for 
weight non-linearity in the beam-forming network. 
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SPECIFICATION 
Adaptive antenna 

5 This invention relates to an accelerated convergence adaptive antenna. 5 
The objective of an adaptive antenna is to select a set of amplitude and phase weights with 
which to combine the outputs from the elements in an array so as to produce a far-field pattern 
that in some sense, optimizes the reception of a desired signal. The substantial improvements 
in system anti-jam performance offered by this form of array processing has meant that it is 
10 now becoming an essential requirement for many military radar, communications and navigation 10 
systems. 

The key components of an adaptive antenna system are illustrated in Fig. 1 . The amplitude 
and phase weights are selected by a beam-pattern controller 10 that continuously updates them 
in respond to the element outputs. The resulting array beam pattern is continuously adjusted to 

15 contain nulls which track the bearings of any interference and jamming sources. 15 
In the illustration, the inputs to the adaptive combiner 1 1 take the form of a primary (or 
reference) channel and a set of auxiliary signals. Usually, some facet of the desired signal is 
utilized to ensure that the primary input contains signals predominantly from this source. In 
addition, the auxiliary channel signals are arranged to provide spatial samples of the jamming and 

20 noise environment only. The multiple channel signals are then combined in order to minimize the 20 
resulting signal power at the beamformed output. Given that the combined signal is expressed 
by 

e(t)-WTx(t)+y(t> (1) 25 

where W is the weight vector and x(y) and y(t) describe the complex baseband envelope of the 
signals at the auxiliary and reference inputs to the combiner, then the residue output power is 
given by: 



25 



30 1 1 

Po = _ E le-ffl e(t)}= — E l(W T x(t)+y(t))* (W^M+ y(t))l (2) 
2 2 



30 



The operator E|.} denotes the expectation function. 
35 The minimum output power condition can be calculated by equating the complex gradient of 35 
this expression to zero. This procedure yields for the solution for W 

Sat-DE|tf«a T «l]- , EU # Wy«l (3) 

40 In this expression, the first factor represents the covariance matrix associated with the signals 40 
received on the auxiliary channels. The second factor is a column vector representing the zero 
time lag correlations between signals received by the auxiliary and reference channels. 

The most commonly employed technique for deriving the adaptive weight vector uses a 
"closed loop" gradient descent algorithm where the weight updates are derived from estimates 

45 of the correlation between the signal in each channel and the summed output of the array. With 45 
the closed loop scheme the summed output from the combiner is used essentially as a feedback 
term to the beam-pattern controller. It is this feature which allows the adaptive control algorithm 
to correct for weight non-linearities and offsets. The closed loop processor can be implemented 
in an analogue fashion using RF correlation loops or digitally in the form of the Widrow LMS 

50 algorithm. ; 

The value of this approach should not be underestimated. Gradient descent algorithms are very 
cost-effective and extremely robust but unfortunately they are not suitable for all applications. 
The major problem with an adaptive beamformer based on a steepest gradient descent process 
is one of poor convergence for a broad dynamic range, signal environment. This constitutes a 

55 fundamental limitation for many modern systems where features such as improved antenna 55 
platform dynamics (in the tactical aircraft environment, for example), sophisticated jamming 
threats and agile waveform structures (as produced by frequency hopped, spread spectrum 
formats) produce a requirement for adaptive systems having rapid convergence and high cancel- 
lation performance. 

60 In recent years, there has been considerable interest in the application of direct solution or 60 
"open loop" techniques to adaptive antenna processing in order to accommodate these increas- 
ing demands. As distinct from the "closed loop" control processor, the open loop techniques 
are termed so because they do not include the sum channel signal as a potential source of 
feedback and hence derive the weight solution directly from the input data samples, in the 

65 context of adaptive antenna processing, these algorithms have the advantage of requiring only 65 
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limited input data to describe accurately the external environment and provide an antenna pattern 
capable of suppressing a wide dynamic range of jamming signals. In fact, these methods usually 
provide an efficient implementation of the control law described by equation (3). The concept of 
rapid convergence rate adaptive arrays is not new and was originally introduced in a classic 
5 paper {"Rapid Convergence Rate in Adaptive Arrays", Reed, I.S., Mallett, J.D. and Brennan, L.E., 5 
IEEE Trans., 1974, AES-10, pp 853-863). Since then many different classes of "open loop" 
adaptive processors have been proposed involving the use of numerical techniques such as Q-R 
decomposition ("Application of a Systolic Array to Adaptive Beamforming", Ward, C.R., Robson, 
A.J., Hargrave, P.J. and McWhirter, J.G., Proc. IEEE, Vol. 131, PE F, No. 6, Oct 1984, pp 
10 638-645), and Gram-Schmidt orthogonalisation ("Introduction to Adaptive Arrays", Monzingo. 10 
R.A. and Miller, T.W., John Wiley & Sons). An approximation to the latter technique known as 
. Sequential Decorrelation has also been considered for this application ("Adaptive Cancellation 
Agreement", McQueen, J.G., UK Patent Specification No. 1-599-035, 30 Sept 1981). 
Open loop algorithms may be explained most concisely by expressing the adaptive process as 
15 a least squares minimisation problem. In fact, the least-squares algorithm defines the optimal 15 
path of adaptation and requires only minimal data to provide a null pattern which cancels 
jamming signals over a wide dynamic range of received power levels. 
The least-squares control algorithm may be derived as follows: 

Consider the simplified block diagram of the adaptive antenna shown by Fig. 1 . The response 
20 from the array at time tj is given by: 20 



(4) 



Here, the complex vector xft) has (N-1) components describing the complex envelope of signals 
25 received across all but one element of the array at time t, while y(tj) describes the corresponding 

signal received at the clamped or reference element. The vector W describes the complex 

weights applied across the array. 
An alternative approach to the least squares estimation problem which is particularly good in 

the numerical sense is that of orthogonal triangularisation, as described by Golub, G.H., "Numeri- 
30 cal Methods for Solving Linear Least-Squares Problems," Num. Math., 1965, 7, pp 206-216. 

This is typified by the method known as QR decomposition which we generalise here to the 

complex case. An n x n unitary matrix Q(n) is generated such that 



35 Q(n)X(n)* 



R(n) 



(5) 



where R(n) is an (N-1) by (N-1) upper triangular matrix. Then, since Q(n) is unitary we have 

(6) 



40 
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' y(n) * 


E{nH|e(n)|H|Q(n)e(n)|H| 




W(n)+ 






0 







where 

45 

u(n)=P(n)y,(n) (7) 
end 

50 y(n)=S(n)y.(n) (8) 



P{n) and S(n) are simply the matrices of dimension (N-1) by n and (n-N+1) by n respectively 
which partition Q(n) in the form 



55 



f D 

Q(n)= — 

L a 



P(n) 
S(n) 



(9) 
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It follows that the least-squares weight vector W(n) must satisfy the equation 



R(n) W(ri>+u(n)=0 
and hence 



(10) 
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Since the matrix R(n) is upper triangular, equation (10) is relatively easy to solve. The weight 
vector W(n) may be derived here quite simply by a process of back-substitution. Equation (10) is 
also much better conditioned since the condition number of R(n) is given by 
5 5 
Cn(R(n))=Cn (Q(n)X(n))=Cn(X(n)) (12) 

This property follows directly from the fact that Q(n) is unitary. 
When the method of Q-R decomposition is applied to signals from an antenna array, the 

10 technique provides the basis for an efficient recursive least squares minimisation process for 10 
adaptive control. When utilized in the "off-line" mode to control a real-time beamformer, how- 
ever, digital beamforming is implicit since the "off-line" processing incorporates no feedback to 
allow for correction of weight non-linearity. 
According to the present invention there is provided an adaptive antenna arrangement includ- 

15 ing a plurality of antenna elements in an array, a beam-forming network, one of the antenna 15 
elements providing a primary signal, the remaining elements providing auxiliary signal inputs to 
the beam-forming network, a beam pattern controller to which the auxiliary signals are applied, 
the controller being adapted to form amplitude and phase weights to be applied to the beam- 
forming network whereby the array beam pattern is continuously adjusted to contain nulls which 

20 track the bearings of unwanted received signals, and means for combining the output of the 20 
beam-forming network with the primary signal, characterised in that the controller derives an 
optimal gradient vector by a least-squares process or an approximation to a least-squares 
process to update the weights which are then applied to the beam-forming network and a 
feedback signal being the array residual signal is applied to the controller to modify the weights 

25 to correct for weight non-linearity in the beam-forming network. A significant feature of the 25 
invention is the aspect concerning the introduction of the feedback signal into the least-squares 
algorithm which enables the system to correct for weight non-linearity and offsets and yet still 
retain a convergence performance which is comparable with an "open loop" control process. 
A preferred embodiment of the invention includes means for time multiplexing the auxiliary 

30 signals whereby the beam pattern controller performs the Q-R decomposition in a timeshared 30 
mode, the auxiliary signals being processed in a cyclic manner. 

Embodiments of the invention will now be described with reference to the drawings, in 
which:- 

Figure 1 depicts key components of a known adaptive antenna arrangement (already referred 

35 to). 35 
Figure 2 depicts the key components of an adaptive antenna arrangement utilising Q-R decom- 
position of the auxiliary signals with feedback of the array residual signal, 
Figure 3 depicts the arrangement of Fig. 2 with reference to a modified Q-R algorithm, 
Figure 4 depicts a modification of the arrangement of Fig. 2 using a timeshared Q-R processor 
40 to control an analogue beamformer, and 40 
Figures 5-8 illustrate the results of computer simulations of the arrangement of Fig. 4. 
By reorganisation of the least squares equations that define the Q-R process, an algorithm can 
be derived incorporating a feedback term. 
At any particular stage during adaptation, the QR or equivalent process provides a computed 
45 weight set W,. to the beamformer 20. Due to non-linearities in the weighting component, this is 45 
translated by the beamformer into a vector Wa such that the output from the network is W a X . 
The computed and applied weight vectors are related by: 



50 



Wc=W,+Aw (13) 

where Aw is some arbitrary weight error term caused by component non-linearity.' However, it is 
known that the least squares processor 21 derives the weight solution Wc by minimising the 
error residual e represented by: 



where 
65 V/-XW.+Y 



50 



55 



55 e-XW^y. (14) 

By substituting the equation for NA^ into (14) we find that: 
e-X(W,+Aw)+y 

60 **-xAw+XW,+y (15) 60 
-XAw+y/ (16) 
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is a measured quantity as indicated by Fig. 2 corresponding to the output from the analogue 
beamformer 20. This equation shows that, by modifying the input data to the QR or equivalent 
process, a weight correction term can possibly be calculated which indicates the error between 
5 the digitally computed weight solution and the weight vector actually applied by the beamformer. 5 
The QR process can be used to minimise the quantity e"e therefore giving a solution to the 
redefined covariance equations: 

(X"X) Aw« -XY (17) 

10 10 
This equation suggests that an iterative procedure could be set up whereby the weight set 
applied to the beamformer, VJk is updated in an update network 22 by the current correction 
term, Aw and the least squares process repeated in order to converge towards a null correction 
(i.e. Aw=0). It would not be necessary for the solution to equation (14) to be computed 

15 explicitly by this method; in fact, W,. would simply be initialised to an arbitrary vector. 15 
Now consider an adaptive antenna algorithm that seeks to minimise the output power from an 
n element array, subject to a fixed unit weight being applied to one of the elements as indicated 
by Fig. 3. If W* denotes the weighting vector applied to the elements with variable weights, the 
expectation of the appropriate steepest descent weight update equation may be expressed in the 

20 form 20 

WK +1 «Wr/iE Irx*} (18) 

Here E irx*( is the correlation between the array output with applied weighting vector W* and the 
25 complex conjugates of the signals at the element outputs with variable weights, the latter being 25 
instantaneously denoted by x. 
If we denote the output from the element with fixed weight by y, we have that 

r=x T W K =Y (19) 

30 30 
Hence 

Wk* ,=Wk-/4E ix # x T l Wk+E |y x'l] (20) 
35 Or 35 
^.i^Wk-^Wk+E |y x"l] (21) 
where 

40 40 
M=E|x*x T { (22) 

is the covariance matrix characterising the element outputs associated with the variable weights. 
The steady state weighting vector, W,,, is the solution of 
45 45 
MV^+Elyx'HO (23) 

We may therefore write the update equation in the form 

50 W^-W.+aMIW.-W/l (24) 50 

or 

[^-^ +1 ]=[^-^//M[^-^] (25) 
55 55 
Let 

AWk+Wo-Wk 

60 denote the error from the steady state weighting vector after the K lh update. We then have that 60 

AWk + t - AWk-atMAWk (26) 

Convergence is now seen to depend on the eigenstructure of the covariance matrix M. If we 
65 modify the update equation by pre-multiplying the correlation E {yx*} by M" 1 such that 65 
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Wk+i^Wk-M -1 E {nc # } (27) 

the corresponding update equation for the weight error takes the form 
5 5 

=[1-/4 AWk (28) 

Convergence is now independent of the eigenstructure of M, and will occur provided |1-/4<1, 
10 corresponding to /z<2 for positive//. 10 
This form of the accelerated convergence update equation is exactly equivalent to that derived 
normally for the more general case of power minimisation subject to an arbitrary linear con- 
straint. The fact that we are dealing with a single fixed weight has enabled the equations to be 
recouched in a form obviating the need to evaluate the pseudoinverse of a matrix. 
15 It is clear from the above equations that 15 

MAWk — E Irx*} (29) 

Hence 

AWk— M~ 1 E |rx*} (30) 

Premultiplying the correlation E {nc*t by to accelerate convergence is thus equivalent to 
calculating the weight error after the K* weight update. 
25 We now consider the method of "QR with feedback" and demonstrate how this technique 25 
produces an extremely effective implementation of the defined accelerated convergence algor- 
ithm. 

Referring to Fig. 3 we have for the summed output from the array, 
30 r-XWi+y (31) 30 

where X represents a block of p data snapshots received across the weighted elements, y 
contains the corresponding p samples of the signal received at the clamped element, Wk is the 
weight set applied at the K lh iteration and r is the vector of the p residual components. 
35 Following the analysis given previously, the weight update, AW* to be computed at the K lh 35 
iteration is obtained by minimising the norm of the residual vector, e where 

e=XAWK+r (32) 

40 The weight update can be obtained by the QR algorithm. For instance, pre-multiplying by the 40 
unitary matrix Q gives 



Q e=Q XAWx+Qr 

45 



in 



AWk+ (33) 



45 



and the least-squares conditon produced when RAWK+b,. 
50 Having thus obtained AWk, the applied weight vector is then updated according to the 50 
equation 

Wkh=Wk+MWk (34) 
55 where is the update gain factor. Substituting for AWk gives: 55 
W^-Wk-zzR 1 b, (35) 

Now, the block estimate of the covariance matrix formed from p signal samples received by the 
60 adaptively weighted elements of the array is given by jft where 60 

M=xhx=rhr (36) 

Assuming jOl, and its associated square root factorisation R, are invertable, then 
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r-i=M-irh (37) . 
Therefore, substitution gives 
5 Wk+i-Wk-^-'R^, (38) 

However, it is noted that 

::ti-[ 

thus giving: 

15 



10 



Qr 
O 



(39) 



20 



However, 



(40) 



R 

Q 



or 



(41) 



25 Therefore, we have 



Wk ♦ i = Wk-/^IvT " 1 X«Q H Q r 



30 



AW 



(42) 



35 



This can be compared directly with the weight update equation describing the accelerated 
convergence technique derived earlier where 



Wkm=Wk~/M iE|x # r} 



(43) 



10 



15 



20 



25 



30 



35 



We notice that M 1 X H r is now simply the block sample estimate of M * E|x # rf. We 

40 therefore conclude that QR with feedback implements precisely this specific form of accelerated 40 
convergence but where we have replaced the expectation operations with block estimates. 

As has been discussed above, for the "Q-R with feedback" system a complete weight update 
was obtained by a full Q-R decomposition of at least N data vectors. Consider now "Q-R with 
feedback" operating in a timeshared mode. Instead of all the signal components present at the 

45 (N-1) adaptively combined channels being fed into the Q-R processor alongside the output from 45 
the analogue beamformer, the array elements are processed in a cyclic manner. A block diagram 
describing this process is shown by Fig. 4. Obviously, decreasing the number of inputs to the 
Q-R processor reduces the system complexity. The performance of such a system-in particular 
the effect on rate of convergence and response to several jammers with a large spread of 

50 power levels-has been assessed by computer simulation. 50 
The simulation was of "timeshared Q-R with feedback" when used in combination with an 
analogue beamformer with a sinusoidal weight coefficient transfer characteristic. Trie program 
allowed cyclic rotation around the (ISM) adaptively combined elements with any number of signal 
components from 1 up to n-1 being taken into the Q-R processor at any one time. 

55 Figs. 5 to 8 demonstrate the adaptive performance for a 5-element array and three indepen- 55 
dent jammers broadly separated in azimuth, of power levels 0, —10, —20 dB relative to a —70 
dB noise floor at each element of the array. 

In Fig. 5 the signal components from each of the adaptively combined elements are fed in turn 
into the Q-R processor, along with the array residue /, to form a weight update term for that 

60 particular element. To ensure sufficient smoothing in the Q-R processor, a dwell time of 10 data 60 
samples on each element was necessary, along with an update gain factor of 0. 1 to ensure 
stable convergence. A time constant of response of —80 samples is achieved against the 
strongest jammer. The rate of adaptation against the weaker jammers is considerably slower. 
In Fig. 6 the size of the Q-R processor has been increased to take in the signal components 

65 from any two elements plus the array residue. The time constant of adaptation against the two 65 
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stronger jammers is now similar (time constant ~60 samples). The response against the —20 
dB damming source shows an improvement over the previous example. 

For the simulation shown in Fig. 7, the size of the Q-R processor now exceeds the number of 
independent jamming sources. The rate of adaptation against the three jammers is identical (time 
5 constant —60 samples). Null depths in excess of 60 dB are achieved against all jammers. 5 
Fig. 8 corresponds to a full Q-R decomposition with alt (N-1) adaptively combined channels 
being fed into the Q-R processor alongside the array residue /. This case is equivalent to "Q-R 
with feedback". Although the rate of adaptation against each jammer is identical, the null depths 
achieved against the jamming signals in the timescale considered, are noticeably less than the 60 
10 dB obtained in Fig. 7. 10 
It can be seen that the system is able to suppress jammers of a wide dynamic range of 
power levels providing the size of the Q-R processor is greater than the number of jammers 
(Ns). If fewer than N s elements are multiplexed, then the rate of adaptation becomes power 
dependent. The technique is particularly interesting since, obviously, convergence performance is 
15 shown to be directly related to processor complexity. For instance, a simple processor (e.g. a 15 
two input Q-R processor) has fast convergence only for the strongest jammer; in fact, the 
processor is essentially performing a power normalised gradient descent algorithm. Increasing the 
size of the Q-R processor allows more and more jammers to be cancelled at equal rate of 
convergence. 

20 It has also been shown that best cancellation performance is obtained when the Q-R proces- 20 
sor size (i.e. the number of input channels) is one greater than the number of jammers. 

CLAIMS 

1 . An adaptive antenna arrangement including a plurality of antenna elements in an array, a 
25 beam-forming network, one of the antenna elements providing a primary signal, the remaining 25 
elements providing auxiliary signal inputs to the beam-forming network, a beam pattern controller 
to which the auxiliary signals are applied, the controller being adapted to form amplitude and 
phase weights to be applied to the beam-forming network whereby the array beam pattern is 
continuously adjusted to contain nulls which track the bearings of unwanted received signals, 
30 and means for combining the output of the beam-forming network with the primary signal, 30 
characterised in that the controller derives an optimal gradient vector by a least-squares process 
or an approximation to a least-squares process to update the weights which are then applied to 
the beam-forming network and a feedback signal being the array residual signal is applied to the 
controller to modify the weights to correct for weight non-linearity in the beam-forming network. 
35 2. An arrangement according to claim 1 wherein the correlation between the summed output 35 
of the array with the applied weighting and the complex conjugates of the signals at the element 
outputs with variable weights is premultiplied by a factor M 1 where M is the covariance matrix 
characterising the element outputs associated with the variable weights. 

3. An arrangement according to claim 1 or 2 including means for time multiplexing the 

40 auxiliary signals whereby the beam pattern controller performs the least-squares process or an 40 
approximation to the least-squares process in a timeshared mode, the auxiliary signals being 
processed in a cyclic manner. 

4. An adaptive antenna arrangement substantially as described with reference to Fig. 2, Fig. 
3 or Fig. 4 of the drawings. 
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